Superfluid Fermi Gases with Large Scattering Length 
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We report quantum Monte Carlo calculations of superfluid Fermi gases with short-range two-body 
attractive interactions with infinite scattering length. The energy of such gases is estimated to be 
(0.44 ± 0.01) times that of the noninteracting gas, and their pairing gap is approximately twice the 
energy per particle. PACS: 03.75.Fi, 05.30. Fk, 21.65.+F 
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In dilute Fermi gases the pair interactions have a range 
much smaller than the interparticle spacing. However, 
when the two-particle scattering length is large, these 
short range interactions can modify the gas properties 
significantly. A well known example is low density neu- 
tron matter which may occur in the inner crust of neu- 
tron stars 0. The two- neutron interaction has a range 
of ~ 2 fm, but the scattering length is large, —18 fm, 
so that even at densities as small as one percent of the 
nuclear density the parameter akp has magnitude much 
larger than one. Bertsch proposed in 1998 that solution 
of the idealized problem of a dilute Fermi gas in the limit 
akp — > — oo could give useful insights into the properties 
of low density neutron gas. 

Cold dilute gases of 6 Li atoms have been produced in 
atom traps. The interaction between these atoms can 
be tuned using a known Feshbach resonance; and the 
estimated value of akp in the recent experiment Q is 
~ —7.4. As the interaction strength is increased beyond 
that for a — — oo, we get bosonic two-fermion bound 
states. In this sense a dilute Fermi gas with large a is in 
between weak coupling BCS superfluid and dilute Bose 
gases with Bose- Einstein condensation j^j. Attempts to 
produce Bose gases in the limit, a/ro — > oo using Fesh- 
bach resonances & re in progress, and their energy 
has been recently estimated using variational methods 
@. 

In the a — > — oo limit k F /m is the only energy scale, 
and the ground state energy of the interacting dilute 
Fermi gas is proportional to the noninteracting Fermi gas 
energy: 



E (p) =£ E FG = £ 
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Baker [?J and Heiselberg [8j have attempted to obtain the 
value of the constant £ from expansions of the Fermi gas 
energy in powers of akp. Heiselberg obtained £ = 0.326, 
while Baker's values are £ = 0.326 and 0.568. 

Fermi gases with attractive pair interaction become su- 
perfluid at low temperature. The BCS expressions in 
terms of the scattering length were given by Leggett , 
and they were used to study the properties of superfluid 
dilute Fermi g function of akp , by Engelbrecht, 

Randeria and Sa de Melo For akp = — oo they 



obtain an upperbound, £ = 0.59, using the BCS wave 
function. These gases are also estimated to have large 
gaps comparable to the ground state energy per particle. 

Here we report studies of Fermi gases with quantum 
Monte Carlo methods using the model potential: 



v(r) 



m cosh (/J,r) 
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The zero energy solution of the two-body Schrodinger 
equation with this potential is tanh(^r)/r and corre- 
sponds to a = — oo. The effective range is 2//i, and in 
order to ensure that the gas is dilute we use pro > 10, 
where ro is the unit radius; pr$ = 3/4ir. All the results 
presented here are for firo = 12; however some of the 
calculations were repeated for pro = 24 and the results 
extrapolated to l/fi — > 0. 

We have carried out fixed node Green's function Monte 
Carlo [HI (FN-GFMC) calculations with trial wave func- 
tions of the form: 



* V (R)=n/( r «')*(R) 



(3) 



where ... and ... label spin up and down parti- 
cles, and the configuration vector R gives the positions of 
all the particles. Only the antiparallel spin pairs are cor- 
related in this \EV with the Jastrow function / (ry)- The 
parallel spin pairs do not feel the short range interaction 
due to Pauli exclusion. 

In FN-GFMC the is evolved in imaginary time with 
the operator e~ Hr while keeping its nodes fixed to avoid 
the fermion sign problem. In the limit r — > oo it yields 
the lowest energy state with the nodes of *£y These 
nodes, and hence the FN-GFMC energies, do not depend 
upon the positive definite Jastrow function. Nevertheless 
it is useful to reduce the variance of the FN-GFMC calcu- 
lation. In the present work we use approximate solutions 
of the two-body Schrodinger equation: 
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v(r) 



f(r <d) = A/(r < d) 



(4) 



with the boundary conditions f(r>d) = l and f'(r = 
d) = 0- The value of d is obtained by minimizing the 
energy calculated using variational Monte Carlo. Note 
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that a dilute Fermi gas is stable even when a — » — oo, 
unlike dilute Bose gases in the a — > oo limit. 

The calculations are carried out in a periodic cubic box 
having pL 3 — N. The single particle states in this box 
are plane waves with momenta k,: 



2tt 



(n ix x + n iy y + n lz z) . 



(5) 



The free-particle energies depend only on / = n 2 , + n 2 + 
n 2 . For N = 14 and 38 we have closed shells having 
states with I < 1 and I < 2 occupied. The commonly 
used Jastrow-Slater (JS) \Jv(R) is obtained by using: 



$S = [A Yl e^}[A Yl e ' ki ' r; ] 



(6) 



in Eq. [21 The more general, Jastrow-BCS \fV(R) has: 

$bcs = A[<j)( ril ,)(f>(r 22/ )...(j){r nn ,)} , (7) 

with n = N/2. The antisymmetrizer A in the $bcs sep- 
arately antisymmetrizes between the spin up and down 
particles. The $bcs describes the component of the BCS 
state with N particles when 
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The nodal surfaces of <&bcs depend upon the pairing 
function 4>(r), and equal those of <&s when Vi = for all 
hi > hp. 

FN-GFMC gives upperbounds to the energy, which 
equal the exact value when the trial wave function has 
the nodal structure of the ground state. Therefore, we 
can determine the <j)(r) variationally by minimizing the 
FN-GFMC energy. We use the parameterization: 

a/exp[ik ; • r] , (10) 

i, I<I C 

0(r) = P(r) + f3(L - r) - 2(5{L/2) , (11) 

0{r) = [l + 7&r] [l-exp[-cbrl] eXP r^ . (12) 

cor 

The function $(r) has a range of L/2, the value of 7 
is chosen such that it has zero slope at the origin, and 
Ic = 4 here. 

The parameters b, c and a/</ c of <fi(r) are optimized 
by choosing a random distribution of initial values and 
measuring the parameters of the lowest-energy (longest- 
lasting) configurations in FN-GFMC calculation. For 38 
particles it produces an optimum set of parameters 6ro = 
0.8, c = 10, a /=0 ,4 = 0.016, 0.466, 0.0068, 0.00091, 0.0007 
which give the smallest FN-GFMC energy having £ = 
0.440(2). Calculations in which f3(r) = give opti- 
mum values a 7=0 ,4 = 0.24,1.0,0.2,0.057,0.035 and \ = 
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FIG. 1: The optimum <j>(r) (solid lines) and the <j>s (dashed 
lines) for iV = 38. Curves ending at L/2, L/\/2 and L-^/3/4 
are in the 001, 011 and 111 directions of the periodic box. 
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FIG. 2: The momentum distribution of particles 



0.459(2), while the Slater </> s (r) having (3{r) = and 
a/=o,4 = 1, 1, 1, 0, gives a much larger £ = 0.54. 

The optimum 0(r) is compared with the </>s( r ) m Fig. 
^ it has a sharper peak at r = 0. This peak depends upon 
the Jastrow function f(r) acting between all the N 2 /A 
antiparallel spin pairs. For example, the 4>(r) obtained 
by solving the BCS equation with the bare potential in 
uniform gas without the f(r) has a much sharper peak. 

The optimum 4>{v) has «o < ol\ for N = 38; in contrast 
the variationally determined BCS wave function has «o > 
U\. The momentum distribution of particles in the trial 
and evolved (r — and 00) wave functions are shown in 
Figure 2. For N = 38 the occupation of the 1 = state 
is smaller than the 1 = 1, calculations with much larger 
values of N are planned to test if this is a finite box size 
effect. 

We have attempted further optimizations by incorpo- 
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rating backflow [12J, [13| into the BCS pair functions <p. 
Initial calculations indicate that this will reduce the £ by 
« 0.02. On the other hand, estimates of the corrections 
due to the finite range of the present interaction indicate 
that going to the — > limit will raise £ by a similar 
amount. Thus our present upperbound for the constant 



£ is 0.44(1). 

In order to estimate the gap A of this superfluid we 
studied differences between energies of systems with odd 
and even number of particles. A general wave function 
with n pairs, u spin up and d spin down unpaired parti- 
cles can be written as: 



$bcs(-R) = A{[cf>{r lv )...(j){r nn ,)] [^i T (r n+ i)...^ uT (r n+u )] [4>n(r {n+1 y)...ip dl (r {n+dy )] } 



(13) 



The unpaired particles are in ip^ and tpji single particle states. We can write this wave function as the determinant 
of an M x M matrix where M = n + U + d |14j . For example, when u = 2 and d = 3 the matrix is given by: 



/ 



0(7-21') 



^il(ri') 
^2i(iv) 
^(ri') 



(r-22') 



<f>(f(n+u)l') 4>{'T(n+u)2' 



fall**) 



n r l(n+d)>) 
0(f*2(n+d)') 



■0ir( r i) 

^lT( r 2) 



■02T(ri) \ 



K r (n+u)(n+d)>) V'lT( r n+«) ^2t(r n +u) 

^i(r( n +i)') 

4>2i(r( n +dy) 

1>n( F (n+d)') 



(14) 



The fact that the general $bcs(R) can be expressed as 
a determinant makes it possible to perform numerical 
calculations for large values of N. When N = 2n, the 
fully paired ground state has u = d = 0, while those of 
systems with N — 2n + 1 have either u or d = 1. 

The FN-GFMC ground state energies for various val- 
ues of N are shown in figure 3. The straight dotted line 
in Fig. 3 is QAAEfg- The calculated energies have the 
odd-even gap expected in superfluids, and well known in 
nuclei. The values of the odd-even gap: 

A(7V = 2n+l) = E(N)-^(E(N-1)+E(N+1)) , (15) 

are shown in Fig. 4. The estimated value of the gap is 
~ 0.9Efg or ~ 2£Efg- In fact the odd particle removal 
energies, E(N = 2n + 1) — E(N = 2n), at fixed density, 
are ~ (4/3)£'fg- The odd particles in the interacting gas 
have energies higher than that for noninteracting. Ap- 
parently the odd particles do not gain any benefit from 
the attractive pair potential, on the other hand they hin- 
der the pairing of the others. BCS calculations including 
polarization correction give A = 0.81Efg in the 

large a limit. 

Several consequences of the strong pairing in this su- 
perfluid gas are seen in the calculated energies. Nonin- 
teracting Fermi gases have shell gaps at N = 14 and 38; 
they are not noticeable in this gas. The ground states 
of 15 and 17 particle systems have momenta with 1 = 1 
rather than the I — 2 in noninteracting states and the 
1 = expected in the limit of strongly-bound pairs. 

Some of the differences between the nodal structures 
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FIG. 3: The E(N) in units of E FG 

of the JS and J-BCS wave functions can be seen by con- 
sidering the case where = Ty . For the JS case, the 
up and down determinants will then be identical and 
the complete wave function will be the square of one of 
these determinants. We now imagine exchanging the po- 
sitions of two pairs by rotating them around their center 
of mass. Since each determinant must change sign, the JS 
wave function must go through zero during this exchange. 
When the pairs are separated by small distances the up 
and down determinants are no longer equal. Thus they 
will change signs at different points along the exchange 
path. We therefore expect a negative region which will 
effectively block these "two-boson" exchanges for fixed 
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FIG. 4: The gap in units of Efg 



as an opening up of the configuration space to allow 
pairs to exchange without crossing a node. For systems 
with a paired ground state, the J-BCS presumably allows 
off diagonal long range order via these pair exchanges. 
In most cases the energy difference between the normal 
state evolved from the JS wave function and the super- 
fluid state evolved from J-BCS is very small (< 0.1%), 
and calculations of the type presented here are difficult. 
However, in dilute Fermi gases with large negative a this 
difference is ~ 20 % and calculations of the supcrfluid 
are possible with bare forces. 

We would like to thank M. Randeria, A. J. Lcggett, 
and D.G. Ravcnhall for useful discussions. The work of 
JC is supported by the US Department of Energy un- 
der contract W-7405-ENG-36, while that of SYC, VRP 
and KES is partly supported by the US National Science 
Foundation via grant PHY 00-98353. 



node calculations. In the J-BCS case, the exchanges can 
occur without crossing a node. In the composite boson 
limit where <p(r) is strongly peaked around the origin, 
there is no sign change under pair exchanges when all 
the pairs are well separated. 

In order to further understand the difference between 
the JS and J-BCS wave functions we studied their nodal 
structure for the following three-pair exchange. In ran- 
domly chosen configurations distributed with \I>j_ BCS (R) 
the three closest pairs ii', jj' and kk' were identified. 
Their center of masses are denoted by S;. The wave 
functions &'(x) are calculated for the positions defined as 
follows: All particles m,m' ^ i,i' , j, j' , k, k' retain their 
positions in the random configuration. The positions of 
i,j, k are given by: 

n = S i + s + x(S j -Si) , (16) 

and cyclic permutations of it. Here s is the relative dis- 
tance between particles in a pair. Those of k' have 
— s in place of s, and the typical value, |s| = 0.25ro is 
used in these studies. The three pairs complete a circu- 
lar exchange ii' — ► jj' — > kk', in the x = to 1 interval. 
We calculate the ratio ^>'(x)/^'(x = 0) for many config- 
urations. Note that *'(l)/*'(0) = 1. In a fixed node cal- 
culation the space where this ratio is negative is blocked 
for the diffusion of the configuration. For JS and J-BCS 
wave functions the ratios are negative, on average, over 
29 and 17 % of the x = to 1 domain. For about half 
of the configurations the J-BCS had positive ratio for all 
values of x, while only 20 % of the JS configurations have 
this property. 

We therefore picture the change in the nodal struc- 
ture in going from the JS to the J-BCS wave functions 
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